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ABSTRACT 

For the Warm Dark Matter (WDM) cosmological model the implications of strongly 
inhomogenous, primordial baryon distribution on sub-galactic scales for Big Bang 
Nucleosynthesis, Cosmic Microwave Background anisotropics and Galaxy Formation 
(including fully non-linear evolution to z=Qi) are discussed, and the inflationary theory 
leading to such distributions is briefly reviewed. It is found that Big Bang Nucleosyn- 
thesis is essentially unaffected relative to SBBN and that the change in recombination 
history at z 1500 — 700 relative to "standard" theory leads to differences in the 
anisotropy and polarization power spectra, which should be detectable by the Planck 
satellite provided systematic eS'ects can be accounted for. Moreover, it is shown by 
fully cosmological, hydro/gravity simulations that the formation of galactic discs is 
only weakly affected by going from smooth to highly non-homogenous, initial baryon 
distributions. In particular, the final disc angular momenta at z=0 are as large as for 
the standard case and the "disc angular momentum problem" is solved to within a 
factor of two or better without invoking (hypothetical) energetic feedback events. A 
very desirable difference relative to the the standard WDM model, however, is that the 
on-set of star (and AGN) formation happens earlier. For the "optimal" free-streaming 
mass scale of Mf ^ 1.5 • 10"'^^ Mq the redshift of formation of the first stars 
increases from z*_i=4-5 to > 6.5, in much better agreement with observational data on 
high-redshift galaxies and QSOs. It will, however, not be possible to push z^,i above 
wlO, because at higher redshifts the gas velocity field is nowhere compressive. Probing 
the "dark ages" will hence enable a direct test of this theory. 

Key words: galaxies: formation — cosmology: theory — cosmology: dark matter — 
cosmic microwave background 



One of the most intriguing challenges in modern cosmology 



ter is elementary particles — note though, that axions, de- 
spite being ultra- light, behave like CDM). Structure forma- 
is determmmg the ongm and properties or the dark matter m . ^-mij • r i r n i i 

° ° i- f- ^jQjj -|-j^g CDM scenario races a number oi well known prob- 

the Universe. Progress can be made either through attempts , , .. i i\ ^ ^ i ii 

, , , , , lems on galactic scales: 1) the steep central cusps problem 

to detect dark matter (elementary particles, massive com- , , , \ „x , , , 

, , , . /, , , .Xtt^ n • n- , , , , n , N (s-g- Moore et al. 1999b), 2) the angular momentum prob- 

pact halo obiects (MACHOs), primordial black holes etc.) ^ ^ n r ^--^ c r, ^- ■ r,r,r,n\ o\ 

J \ /I r- / jgjjj (g g >3ommer-Larsen, Gotz & Portman 2003), 3) the 

or indirectly by studying the implications of various dark . . , , , / t^i ■ ^ i ir,r,r,\ i 

, , r 1 r ri r missiug satcllites problem (e.g. Klypm et al. 1999) and pos- 

matter candidates for the formation oi galaxies, clusters oi ., i ^\ i. i .. i i / . i -.^^^ \ 

, sibly 4) the disc heating problem (e.g. Moore et al. 1999a). 
galaxies and the large-scale structure. The former approach ■ j_ j_i i i r-itTT^T,i\ j_ j_ r 

° ■ J c • c -u Going to the warm dark matter (WDM) structure formation 



has so far not given any definite results (except for the pos- 
sible detection of micro-lensing by ~O.5M0 MACHOs in the 



scenario can considerably improve on problems 2-4 and pos- 

^ , . , , , . , , r . n , , sibly also on problem 1. On much larger scales, the scarcity 

Galactic halo), so one is currently lett with the latter. „ , . . . , , , , j_ j_- ^ i i r m-Mvj- 

' or galaxies m voids could be a potential problem for CDM 

The salient feature about conventionar cold dark mat- r,ooi\ i - i i u j- j u • ^ 

T\ , , , , n n iPeeble&.2001i). which also may be r emedied by going to 

ter (CDM) is that as long as the dark matter particles are WDM iBode. Ostriker fc Turokjlifflih . 
much heavier than 1 keV, the actual particle mass does not 

The WDM structure formation scenario has, however, 
also at least one major problem in that it is difficult to 
get early enough galaxy formation to match observational 
* E-mail: jslarsen@tac.dk constraints from high-z galaxies and QSOs. For a free- 
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strea ming: mass scale of Mf ^ 1.5 10^^ h^^ Mq, 
which ISommer-Larsen fc Dolgovl (1200111 found was the op- 
timal value for solving the angular momentum problem, the 
first star formation takes place at 2:<6 (see section 5). In 
contrast, both galaxies and QSOs have now been observa- 
tionally detected at z > 6 (e.g., Hu et al. 2002, Fan et 
al. 2001) and the recent WMAP results possibly indicate 
that the rcionization of the Universe took place at z^c ^ 10 
(IS^^rgslot al. 200:^. 

The implicit assumption in standard WDM galaxy for- 
mation simulations is that at early times the baryons and the 
dark matter are distributed in the same way. The possibility 
exists, however, that the baryons could have a strongly in- 
homogeneous initial (primordial) density distribution. If on 
the (small) scale of these "baryon" clouds the fluctuations 
in the total matter density (baryons, dark matter and ra- 
diation, including neutrinoes) are isocurvature fluctuations 
then the baryon clouds will not be smoothed out by Silk 
damping, nor will they suffer growth by gravitational insta- 
bility a ny faster than the fluctuat ions in the dark matter 
density ijNaselskv fc Novikovll20o3) . Because of the higher 
gas density in the baryon clouds, the clouds will however 
cool out and form stars earlier than is the case in standard 
WDM simulations. The purpose of this paper is to quantify 
the effects of such initially inhomogeneous baryon distribu- 
tions on the epoch of first star/galaxy formation as well as 
the subsequent disc formation. We also calculate how Big 
Bang Nucleosynthesis (BBN) and Cosmic Microwave Radi- 
ation Background (CMB) anisotropics are afi^ected by going 
to this scenario, which we shall denote the WDMB model 
in the following. 

In section 2 we describe the mix between adiabatic and 
isocurvature perturbations in the WDMB model, in section 3 
inflationary models leading to such initial baryon clouds and 
in section 4 the resulting BBN and CMB anisotropy and 
polarization for the WDMB model. In section 5 we present 
our fully cosmological, hydro/gravity disc galaxy formation 
simulations, and finally, in section 6, we summarize our con- 
clusions. 



2 THE MIX BETWEEN ADIABATIC AND 

ISOCURVATURE PERTURBATIONS IN THE 
WDMB MODEL 

In the framework of modern theories of inflation there 
are a lot of possibilities for having a mix of perturba- 
tio n modes in the composit e fluid, which are discussed 
bvlRiazue lo fc Langloi^ ('200rf), 'Bartolo. Matarrese & Kioto' 
|200lh , JPolarski fc Starobinskv (1994), Abramo & Finclh 
J200lli . ISucher. Moodlev fc Turokl 1I2OOIII and others. The 
general idea for classifying modes of perturbations is based 
on a simple deflnition of the isocurvature modes. They do 
not perturb the gravitational potential meaning that the 
initial fluctuations of the total matter density ptot are zero. 
In terms of the total density perturbation Sptot this corre- 
sponds to 



5pu 



N 



where pi denotes the density of each massive species includ- 
ing baryons and different kinds of CDM and WDM parti- 
cles, 5i = 5 Pi/ Pi is the density contrast for each massive 
component, i?^^ is the density ratio between neutrinos pi, 
and black body radiation p-,, and ST /T are the CMB tem- 
perature perturbations. 

Using Eq. Q one can flnd a certain peculiar mode 
(or modes) which compensates the contribution of bary- 
onic perturbations to (5ptot, i.e. it corresponds to the con- 
dition pb^b = —pdmSdm (we assume here without loss of 
generality that there is just one main dark matter compo- 
nent). We will denote this mode a co mpensating isocurva- 
ture mode (CIM) for the dark matter iNaselskv fc Novikovl 
200^. As one can see this mode corresponds to 5T/T = 0. 
This means that the CIM are equivalent to an isotempera- 
ture pertu rbation. Note that exactly the same mode was de- 
scribed bv I Abramo fc Finelli i200l[) . but for compensation 
between quintessence scalar fleld perturbations and some of 
the CDM particles. 

The new aspect in the evolution of the isocurvature per- 
turbations comes from the WDM model. If the rest mass of 
the dark matter particle is comparable to rrix — 1 keV, then 
practically all perturbations at spatial scales smaller than 
the free streaming scale Afe ~ 0.2 ( gwDM.dcc -^-i/s (^^y^ 
Mpc are damped while at scales above Afs they behave like 
in the CDM model. The corresponding froo-streaming mass 
in terms of the scale Afs is fSommcr-Larscn fc Dolgov 200:0) 



Mfe = 3.7 10" fiwDM h'^ 



0.1/i-iMpc 



Ma 



(2) 



.ST 



^ p^S^ + 4p-,{l + R„^)— =0, 

i=0 



(1) 



where h = _H'o/(100km/s/Mpc) is the Hubble constant, 
gwDM,dec is the effective number of particle degrees of free- 
dom when the WDM particle (s) decouple from the cosmic 
plasma, and 

n^uMh = ( ) (j^j • (3) 

This means that for pure adiabatic initial perturba- 
tions in the WDM model all inhomogeneities at the scales 
M < Mfs are damped and baryons and WDM particles fol- 
low an approximately homogeneous and isotropic spatial dis- 
tribution. 

The situation changes dramatically if we instead of pure 
adiabatic perturbations have a mix between adiabatic and 
isocurvature modes of fluctuation. Let us discuss a model 
with isocurvature baryonic perturbations which have a typ- 
ical spatial scale A ^ Afs, but the corresponding amplitudes 
of perturbation are much higher than for the adiabatic tail 
of perturbations. This model is close to the baryon isocurva- 
ture fluctuations model (Dolgov fc Silk 1993) and manifests 
itself as a combination of an adiabatic WDM model plus 
small-scale baryonic clouds. We shall denote this model the 
WDMB model. 

Although we will consider in the non-linear galaxy for- 
mation calculations (Section 5) a strongly in-homogeneous, 
initial baryon density fleld -|- CIM, within the framework 
of isocurvature modes even if the the condition phSh = 
—pdmSdm is not satisflcd (e.g., with initially strongly non- 
homogenous baryon -I- smooth dark matter density fields) 
the baryon density fleld will neither be smoothed out due to 
Silk damping, nor grow non-linear much before the dark 
matter density fleld due to gravitational instability (e.g. 
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Naselsky & Novikov 2002). We note, however, that for 
baryon "cloud" masses Md < Mdut — 10 — 20 Mq the baryon 
density field will be smoothed before the epoch of recomb i- 
nation due to baryon dif fusion iZeldovich fc Novikovlll983l) . 

Finally we note that lDemianski fc Dorosh kevichI ll2003ll 
recently also argued for small scale peculiarities in the power 
spectrum, albeit on other grounds. 



3 BARYONIC CLOUDS AS A RELIC OF THE 
BARYOGENESIS EPOCH 

It is necessary to note that the idea of a non-homogeneous 
distribution of baryonic matter at small scales is not new. 
The importance of entropy perturbations in the history 
of the cosmological expansion was ad hoc de monstrated 
by iDoroshkevich. Zel'dovich fc Novikovl a^nd Pee- 

bles (1967,1994) and recently generalized taking into 



the multi-species structure of the cosmological 


by Gncdin & Ostriker 


l|l992f). 'Hoffan & Loebl 


Peebles & Juszkicwich 


(299^. The possible inho- 



plasma 
(1993), 

mogeneities in the baryon distribution at the epoch of 
nucleosynthesis (Inhomogeneous Big Bang Nucleosynthesis 
(IBBS)) has been widely discussed in the literature (see 
the review by Jedamzik & Rehm 2001) in connection with 
the quark-hadron phase transition. But the typical scales 
of such inhomogeneities are extremely small compared 
to the typical mass scales of interest here - see below. 
Other events or processes have been suggested as possible 
sources of isocurvature perturbations partly connected 
with baryon re-distribution in space. For example, cosmic 
strings and corresponding currents and magnetic fields 
could generate s pecific features in baryonic matter distr i- 
bution j Malancv fc Butlei|ll989tl.lYokovama fc Satol l|l99l " 
^ Pohtfsk^^^^ robingk^ "^ 



DoIeov & Silk ( 199^ 
Novikov. Schmalzing fc Mukhanov 



11994 

have suggested 



various ways of generating isocurvature perturbations in the 
frame work of inflation theo ry. Re c entlv l Naselskv fc Novikovl 
l|2002h and pDoroshkcvicli et alJ i2003h have shown that 
baryonic clouds could mimic the acceleration of hydrogen 
recombination at redshift z ~ 10'^ and produce corre- 
sponding peculiarities in the cosmic microwave background 
(CMB) anisotropy and polarization power spectrum. 

Probably, the most interesting question in modern 
baryogenesis theory is how exactly the primordial soup 
of particle species became charge symmetry breaking and 
what the spatial distribution of the baryons might be. We 
will consider one of the most attractive models for bary- 
on ic cloud creatio n just after the end of infla tion, suggested 
bv iDolgov fc Silkl Jl993l) and similar to the lAffeck fc pTii^ 
l|l985fl barvogenesis model. We will assume that two scalar 
fields $ and play as crucial a role in inflation as in charge 
symmetry breaking. Namely, that the field $ determines the 
general properties of inflation while spontaneous CP viola- 
tion is achie ved by a condensate o f the complex scalar field 
(j). Following iDolgov fc Silkl lll993f) . we will describe charge 
symmetry breaking using as the potential of the (j) field 



V(0) 



where 

2 

m^B = 



-ml + PT^ 



/i($ - $c 



(4) 



(5) 



is the effective mass of the (j> field, including interaction with 
the inflaton $. mo is the vacuum mass of 0, "I>cr the critical 
value of the $ field when the effective mass rn^g has a global 
minimum, and the PT^ term corresponds to the temperature 
correction arising from the interaction of the (j) field with 
the thermal bath, assuming no conformal coupling between 
(j) and the curvature R. That all the constants mo. A, /i, <l?cr 
and /3 are assumed to be free parameters of the model allows 
us to optimize the resulting baryonic cloud properties. 

Let us briefiy describe some important stages in the $ 
and (/!> field evolution during and right after the end of infia- 
tion. When inflation started, the initial value of the $ fleld 
is assumed to be higher than the Planck mass, $in > mpi, 
and higher than $cr as well, and m^g > 0. The poten- 
tial V{(l>) has a global minimum ai (p = and roughly 
V{4>) oc m^g|(^p. As one can see, the charge symmetry is 
unbroken. During the slow roll epoch of inflation the ampli- 
tude of the <& field is a decreasing function of time, which 
is typical, for example, for the model of chaotic infiation, 
and when $ — <I>cr the effective mass of the (p field reaches 
the minimum value m^g — —m^ + PT^ < 0. The end of 
infiation corresponds to the condition $ ~ <C ^cr, while 
m^ff = —mo -|- PT^ + fJ-^cr > 0. Obviously, in the vicinity of 
the point $ = the effective mass m^g ^^j^, can be nega- 
tive, or positive but small, which means that it may happen 
that for some period of time |m(,fi,min| < Hst, where -ffsr 
is the value of the Hubble parameter durin g the slow roll 
decrease of the <E> field caused by inflation. iDolgov fc Silkl 
(1993) have shown that the phase transition of the (j> field 
during the epoch |meff,min| < Hsr is essential for the forma- 
tion of bubbles, which are governed by the quantum fluctu- 
ation of the (j) field. The corresponding nu mber density of 
the b ubbles as a function of their mass M is ijDolgov fc Silkl 
Il993h: 



dn(M) 
dM 

where 



M^ 



exp 



7 , 2fM\ 
Mo — ^rnpiH'^ exp [2J/'sr(fc — 



(6) 



87" 



and a — S + I67 



5 = 



7 = 



47r^/i /d<I>y 



tc means the end of infiation and ta is the moment when 
^ = 'I'cr. For the $ field using slow roll approximation we 
obtain 

1 dVinm 



d$ 



(7) 



where Vin{^) is the infiaton potential. That lead to the fol- 
lowing value for the 7 parameter: 



7 = 



47rV f dV^„($) 



27XHg \ d$ 



(8) 



As one can see from Eqs. JSJ - ijHJ, the distribution 
function dn{M) /dM is strongly peaked at M ~ Mo, if the 
corresponding value of the 7 parameter is 7 3> 1. In this 
case we can approximate dn(M)/dM oc S{AI — Mo), where 
S is the Dirac (5-function, and the corresponding mass Mo 
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is a functi on of the fundamenta l parameters, namely, Hst, 
te and tcr- iDolgov fc Silld (Il993l) have shown that the sim- 
plest way to have a large value for the baryon asymmetry 
inside the bubbles is to assume that CP-nonconservation is 
proportional to the value of the scalar field which is much 
larger than that outside the clouds. Thus, because of bub- 
ble formation during inflation, and using CP-nonconserving 
effects, the baryonic matter can have some strong perturba- 
tion strongly peaked at some characteristic mass scale Mq. 
It is not possible at present to calculate the value of Mo 
from first principles, but we shall assume in the following 
that Afdiff < Mo < Mfe (Section 2). 

For simplicity we suppose that all baryonic clouds have 

1/3 

the same characteristic size _Rci oc Mq , which is much 
smaller than the size of the horizon 7?rec at the epoch of 
recombination {z ~ 10^), Rd ^ -Rroc- We denote by pb.in 
and pb,out the baryon density inside and outside the clouds. 
Then the mean density pb,mcan at scales much greater than 
_Rci and distances between them, is 

Pb,mcan = Pb,in/ + Pb,out(l ^ /) , (9) 

where / is the volume fraction of the clouds. We denote 

Pb.iii 



V 



pb.out 



(10) 



Using this one obtains the following relations between 
the mean value of the baryon density and its inner and outer 
values 



Pb.in 



and 



Pb,out 



VPb,, 



1 + /('?-!) 



Pb,mcan 

1 + fiv - 1) ■ 



(11) 



(12) 



Using the quantities / and r; we can define the baryonic mass 
fraction of the clouds 



Vf 



i + fiv-i)' 



(13) 



Obviously, all the parameters /, and Fh are the result 
of fine tuning the inflaton Vin{^) and (j) scalar field potentials 
leading to the formation of the baryonic asymmetry in the 
Universe. 



4 BARYONIC CLOUDS AND THEIR 
COSMOLOGICAL CONSEQUENCES 

There are several sources of information on Sib ~ ph/pa, 
where pb and pcr are the present values of the baryonic and 
critical densities. One comes from the confrontation of the 
Standard Big Bang Nucleosynthesis (SBBN) theory with ob- 
servational data (see the review by Fukugita, Hogan & Pee- 
bles 1998) . The corresponding value of the baryonic density 
from this line of argument is fib/i^ = 0.019 ± 0.001. An ad- 
ditional empirical relation between baryonic and total mat- 
ter fractions fh = Qb/^M ~ 0.1 — 0. 15 for h = 0.7 comes 
from X-ray data on clu sters of galaxies dCarlberg et al.ll996l : 
lEttori &c Fabianll999t) . For the most popular ACDM cosmo- 
logical model with S7m — 0.3 and SIa — 0.7, the correspond- 
ing value of Q^h^ is ^ 0.02, in agreement with the SBBN 
predictions jMason et al.ll200^ . 

An independent and important information about the 



baryonic fraction in the Universe come s from recent CMB 
exper iments such as BOOMERANG ide Bernardis et al 



20(30|) . MAXIMA- 1 jHananv et alj|2 000h. CBI jMason eral 



12003), and DASI jLeitch et alJ 120021. Fitting the CMB 
anisotropy power spectrum to the above-mentioned obser- 
vational data (Tegmark fc Zalda rriaga 2000; Whit e et alJ 
l2000l : lL'esgourgues fc Pelosi^ll20o3) indicates that the baryon 
fraction parameter should be significantly larger than 
the SBB N expected value, namely, Q^h^ — 0.03. 
However, iBond fc CrittendenI (1200111 showed that new 
BOOMERANG, MAXIMA-1 and DASI data do not con- 
tradict fibfe^ = 0.022 ± 0.004. 

It is worth noting that the above mentioned methods 
for the baryonic fraction density estimation from CMB and 
SBBN predictions are based on the simple idea that the 
distribution of matter (including dark matter particles and 
baryons) is nearly homogeneous on all scales, except for 
small fluctuations leading to galaxy and large-scale structure 
formation. Typically, these are assumed to be adiabatic. One 
can ask how sensitive conclusions based on the BBN -I- CMB 
data are to the presence of small-scale baryonic (non-linear) 
clouds before cosmological recombination. 

4.1 Big Bang Nucleosynthesis (BBN) 

In the framework of the SBBN model the baryonic frac- 
tion of matter and its spatial distribution play a crucial role 
starting from the epoch, when the balance between neu- 
trinos (vcT^), neutrons (n) and protons (p) in reactions 
n + Uc ^ p+e~ , n + e'^ ^ P+T^i n p + e^ was broken. 
The corresponding time of violation of the neutrino-baryon 
equilibrium is close to T^^,p — 1 sec when the temperature 
of the plasma was close to Ti,^^p ~ 10^" A" (see the review 
by Olive, Steigman & Walker 2000). The time scale Ty„,p 
determines the characteristic length Z^^.p ~ cr^^.p, which in 
terms of the baryonic mass fraction of matter corresponds 
to 



mpi 



310" 



0.02 



Mo, (14) 



where tpi is the Planck time, pb and p-, are the densities of 
baryons and radiation in the standard cosmological model 
without baryonic clouds. The light elements (chiefly, ''He 
and deuterium) were completely synthesized by Tend ~ 3 • 
lO'^ — 10^ s. In terms of the baryonic mass scale it corresponds 
to 

«„..«.,(^)"Lie(S)"Y%gV%^(«) 



0.02 



Thus, if the characteristic mass scale Mo for the baryonic 
clouds is higher than Mend (~ Mdifi cf. section 2), the cos- 
mological nucleosynthesis in each cloud is independent of the 
others and the mean mass fraction of each chemical element 
is related with the mass fraction of the clouds as 

V J \ pmoan / 



Pout(r) 

pmcan 



Vft""' (pout(r)) 



(16) 



where Vd/V = f, V is the volume and Y^'^, Y^""* are the 
inner and outer mass fractions of the fcth element. If we 



© 2003 RAS, MNRAS 000.1110 



Y P 

r 



High-z Galaxy Formation with WDM 5 



Y „ ( X 10 




10 1 10 



4 



2 



1 

V 

Y P 




1 10 
V m 



Figure 1. The mean mass fraction of *He as a function of »7ni and 
rj for different baryonic mass fractions _F^,. From top to bottom 
Fb = 0.2,0.5,0.8. 

assume that inside the clouds the density pin{r) is uniformly 
and isotropically distributed then from Eq. II1H - 113II and 
Eq. 1161 we obtain 

(n) = n"(pin)Fb + Yr\pout){i - Fb) (17) 

In Fig. and Fig. |5| we show the dependence of the mean 
fraction of "He and D on 77^ = 10^°(nb)/n^ = 274{r2b)^^ 
and 77 for such a model. The three panels of Fig. Q show 
different values of the mass fraction Fh = 0.2 (top), -Fb = 0.5 
(middle) and -Ft, = 0.8(bottom). As one can see from Fig.0 
if the {fth)h'^ parameter is close to (fib)'*^ ~ 0.022 (rim ~ 5), 
and rj = 1, then the mean mass fraction of helium is equal to 
Yp in the standard model (without clouds) and corresponds 
to Yp = 0.247. For the model = 0.2, 1 < < 10 and 77^ ~ 
5 the mean helium mass fraction lies in the range 0.247 < 



V m 

Y B ( X 10 




Figure 2. The mean mass fraction of D as a function of rjm and 
T) and the same _Fb parameters as in Fig. 1. 



(Yp) < 0.252 which corresponds to the 2% relative deviation 
from the value Yp = 0.247. For the model with Fb = 0.5, as 
it follows from Fig. the maximum value of (Yp) = 0.256 
is achieved for = 10 and the corresponding deviation is 
3%. For the model with F^ — 0.8 we obtain a quite similar 
result: (Yp) < 0.256 for all values in < 10 . Since the 
uncertainty in the observationally determined value of Yp 
is ~ 5% (when systematics are included — B. Pagel 2003, 
private communication) then even if rjm was known to high 
precision (e.g. from the CMB data) it would not be possible 
to distinguish between the smooth and clumpy initial baryon 
distribution scenarios using BEN. 

For the primordial mass fraction of deuterium the de- 
pendence of (Yo) on and -Fb is shown in Fig. |21 As 
can be seen from the figure, for rjm ~ 5 and ri — 1 the mass 
fraction of deuterium is close to (yo.st) — 4- 10~^. Moreover, 
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for 77 < 10 and Fb = 0.2, 0.5, 0.8 the corresponding variation 
of (Id) is < 20 — 30% for rjm ~ 5. As this is less than 
the uncertainty on the observationally determined value of 
(^0) (~ 40% - B. Pagel 2003, private communication) this 
can not be used either to select between the two classes of 
models. 

For both "He and D, if 77 is increased beyond 77=10 the 
further change in Yp and Yo for a given 77ni is small compared 
to the above mentioned observational uncertainties. 



4.2 CMB Anisotropy and Polarization 

As was mentioned by 'Nasel sky fc Novikovl ||20o2), the pres- 
ence of the baryonic clouds in the primordial hydrogen- 
helium plasma at redshift z ~ 10'^ changes the dynamics of 
hydrogen recombination due to the non-linear dependence 
of the ionization fraction Xc on the baryon density. During 
the period of recombination diffusion of baryons from inner 
to outer regions of the clouds can suppress any small scale 
irregularities inside the clouds. The natural length scale of 
this process is of the order of the Jeans length Rj ^ CsUec, 
where Cs is the baryonic sound speed and tree is the corre- 
sponding time when the plasma bec ame transparent for th e 
CMB radiation kUu et al..,200L : .Naselsky fc NovikovlbOOSl) . 
For adiabatic perturbations at scales M ^ 10^^ M© (which 
are the sources of the Doppler peaks in the CMB anisotropy 
and polarization power spectra) the evolution during the pe- 
riod of recombination depends not on the ionization fraction 
inside or outside the clouds, but rather on the mean value 
of the ionization fraction at the scales of adiabatic pertur- 
bations. This mean ionization fraction does not correspond 
to the ionization fraction for the mean value of the baryonic 
density due to non-linear effects. 

For the cloudy baryonic model we introduce the mean 
value of the electron density 



(ne) = Tie, in/ + 7Jc,out(l — /), 



(18) 



where nc,in and nc,out are the number densities of free elec- 
trons inside and outside the clouds and 



^c,out 



(19) 



where x^in, a^o.out, yncin and yncout are the ionization frac- 
tions and helium mass fractions for the inner and outer re- 
gions, and the nucleon number density is denoted Uh- Note 
that by definition Xc = tIc/tih, where tih is the number den- 
sity of neutral and ionized hydrogen. Let us introduce the 
mean value of the ionization fraction (xe), 



{:r.) = M(l-(yHo))-\ 
then 

where 

JJ ( 1 - Yp,, 



(20) 



(21) 



Gin — 
Gout 



1 + /(^ - 1) V 1 - (Yp) 

1 — / I 1 ~ Yp^out 



and (Yp) denotes the mean mass fraction of helium. 



(22) 
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Figure 3. The rate of ionization for the clumpy baryonic model. 
The dash-dotted line corresponds to the ionization fraction for 
the inner regions. The dotted line is the ionization fraction for the 
outer regions. The solid line is the mean fraction of ionization. The 
dashed Hne is the ionization fraction for the model with fi^/i^ = 
0.022 without clouds. 




In Fig. 13 we plot the functions Xc,in, Xa 



and (a;e) 



for (fib/i^) = 0.022 , h = 0.7, 77 = 11, and / = 0.1, and 
cosmological parameters SIk = (the curvature parameter), 
r^A = 0.7 (vacuum density) and riwDM/i^ = 0.125 (WDM 
density) . A modified version of the RECFAST code was used 
for the calculations. 

As one can see from this figure the fractions of ioniza- 
tion (xc), Ko.in and Xc.out have different shapes, and differ- 
ent asymptotic behaviours at "low redshifts". In the range 
700 ^z< 1500 recombination is somewhat accelerated and 
for z < 700 slightly delayed compared to "standard" recom- 
bination. 

In order to compare the CMB anisotropy and polar- 
ization power spectrum in the cloudy baryonic model with 
the st andard one (without clouds ) we modify the CMBFAST 
code iSeliak fc Zaldarriagall993) . by taking into account the 
more complicated ionization history of the plasma. For the 
comparison we assume in both cases (late) full reionization 
corresponding to a value of the optical depth of — 0.1. 
is related to the redshift of full reionization Zr as 

-2/3 



(l-{Yp) 



0.022 



-2/3 



0.76 



0.125 



1/3 



(23) 



so Tr — 0.1 corresponds to a redshift of reionization of 
Zr ~ 14 (the following results do not change in any sig- 
nificant way if rr=0.05 or 0.15 is adopted instead). In Fig.0] 
we plot the CMB anisotropy power spectra for the stan- 
dard and cloudy WDM models together with the latest 
BOOMERANG, MAXIMA-1, and CBI observational data 
(the very recent WMAP data are essentially in complete 
agreement with the above data, Spergel et al. 2003). As can 
be seen from the figure small deviations of the power spec- 
trum appear, but both models are in agreement with the 
data of the CMB experiments (see also Table 1). 

To illustrate the effect of baryonic clouds in the WDM 
model and to compare with the standard WDM model, tak- 
ing into account the expected sensitivity of the PLANCK 
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Figure 4. The CMB power spectra for the WDM models. The 
sohd hne corresponds to the standard model, the dashed line to 
the cloudy baryonic WDM model. 
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Figure 5. The CMB polarization power spectra for the WDM 
models. The solid line corresponds to the standard model, the 
dashed line to the cloudy baryonic WDM model. 



experiment, we describe the deviations of the anisotropy and 
polarization power spectra in terms of the functions 

D^l) = 2 [C^{1) - C?(/)] / [CS{1) + C(0] , 

D^(l) = 2 [C^{1) - C^l)] I {CliV) + Cl{l)\ , (24) 

where (l) and C" (1) are the CMB anisotropy power spec- 
tra for the baryonic clouds and smooth baryon models, re- 
spectively, and (Z) and Cf(Z) are the corresponding polar- 
ization power spectra. Obviously, it is necessary to compare 
the functions D°'{1) and D''{1) with the errors of the C{1) 
extraction for the PLANCK mission. 

Assuming that systematic effects can be essentially re- 
moved, the uncertainty should be close to the cosmic vari- 
ance limit 
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Figure 6. D°'{1) and D''(/)for the CMB anisotropy and polariza- 
tion power spectra in WDM models. The solid line corresponds 
to anisotropy, the dash line is the polarization. The thick solid 
line corresponds to the cosmic variance effect. 



Wi ~ exp 



2li 



(25) 

/sky — 0.65 . 



Here /sky is the sky coverage during the first year of 
observations, cTp is the sensitivity per resolution element 
^'fwhm X 6'fwhm and Is = a/8 In 26'p^jj^ . In Fig. |Slwe plot 
the functions £)"(/) and D^{1) for 2 < Z < 2000 and show 
the uncertainty due to cosmic variance. 

For example, for the HFI PLANCK 214 GHz chan- 
nel ^FWHM is 5.5 arcmin and the errors at I < 2000 are 
dominated by the first term in Eq. 1251 . As one can see 
from Fig. |S| at Z ~ 1000 — 1500 the peak-to-peak ampU- 
tudes are ~ 10 - 20%, while AC{l)/C{l) ~ 5 - 7%. This 
means that such signals in the anisotropy power spectrum 
can be detected by the PLANCK mission, if the systematic 
effects can be accounted for to better than ~ AC{l)/C{l) ~ 
[/sfe,(Z + l/2)]-i/^ 

As one can see from Fig. |S|the polarization power spec- 
tra are even more sensitive to distortions of the ionization 
history of the Universe in the WDM mode l with bary- 
onic clouds. Recently the DASI experiment iLeitch et alJ 
2OO4) has detected CMB polarization without a significant 
presence of polarized foregrounds. This implies that the 
PLANCK mission potentially will be able to measure the 
polarization of the CMB with unprecedented accuracy, prob- 
ably close to the cosmic variance limit, which means that 
peculiarities of the polarization power spectrum caused by 
baryonic clouds could be observationally detected, if present, 
during the coming years. 



5 DISC GALAXY FORMATION WITH 
BARYONIC CLUMPS 

The uneven distribution of baryons in this model could of 
course have consequences for the formation of disc galaxies. 
The baryonic clumps are additional structures, which would 
not be present in "normal" galaxy formation scenarios. They 
could possibly disturb the assembly of discs, similar to the 
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Table 1. CMB power spectrum parameters for the WDM and WDMB models. 



Obs/Mod 


BOOM 


MAXIMA 


CBIMl 


CBIM2 


VSA 


Dof 


12 


10 


7 


7 


20 


Standard model (WDM) 


11.2 


14.1 


3.20 


7.15 


7.08 


Cloudy model (WDMB) 


8.69 


11.3 


1.60 


5.62 


6.59 



angular momentum problem in CDM cosmologies (see, e.g., 
Navarro & Benz 1991, Navarro & White 1994, Navarro & 
Steinmetz 1997, Bryan 2002, Sommer-Larsen, Gotz & Porti- 
nari 2003). For CDM, the early formation of small haloes 
leads to many merging events, during which the subclumps, 
which later form the central disc, can loose energy and or- 
bital angular momentum by dynamical friction, resulting in 
discs with too little specific angular momentum. The pres- 
ence of baryonic clumps may aggravate this problem, and 
may even effectively prevent the formation of galactic discs. 
To investigate this, we ran a series of TreeSPH simulations of 
galaxy formation, described in this section, with and with- 
out such clumps. 

5.1 The Code 

We use the gridless Lagrangian N-body and Smoothed 
Particle Hydrodynamics code Tree SPH described in 
ISommer-Larsen. Gelato fc Vede _^^99|), which is modelled 
after that of lHernauist fc Katd il989l') . 

Radiative gas cooling is included in the simulations 
JSutherland fc Dopitalfl993h . Some are run with a primor- 
dial abundance cooling function, in which case radiative 
heating due to a metagalactic UV field is included. This 
homogeneous and isotropic, redshift-dependent UV back- 
ground corresponds to the UV field pr oduced by AGNs 
and y oung galaxies. It is modelled after iHaardt fc Madaul 
il996h and switches on at a redshift of z = 6. Its ef- 
fe cts on the cooling function are included a s discussed e.g. 
in lVedel. Hellsten fc Sommer-LarsenI lll994l) . Other simula- 
tions were run with an about 1/3 solar metallicity cooling 
function (more precisely with [Fe/H] = —0.5). This is the 
metal abundance of the intracluster medium and can prob- 
ably be considered a reasonable upper limit to the metal 
abundance of disc forming gas. In this case, the metagalactic 
UV field is not included, since it does not play a major role 
because of the very high cooling efficiency due to the metals. 
Furthermore, all simulations incorporate inverse Compton 
cooling, which is also explicitly redshift-dependent. 

For the gas dynamics, the smoothing length of each SPH 
particle is adjusted to keep the number of neighbors close 
to 50, and the shear-free Balsara viscosity (Balsara 1995) is 
used, instead of the standar d Monaghan-Gingold viscosity 
jMonaghan fc GingoldllTgsll . 

5.2 The Star Formation Recipe 

Star formation is accounted for by converting SPH particles 
into star particles, keeping the total number of these parti- 
cles constant. We assume that star formation sets in when 
the hydrogen number density exceeds a threshold of nu,c~ 
0.01 cm~^. This may seem quite low, but is still large enough 



to ensure that the gas has cooled below a critical tempera- 
ture of T ~ lO** K, at which the radiative cooling function 
is effectively truncated. In addition to the density criterion, 
the divergence of the velocity field has to be negative in 
the star forming region. As will be seen later, this second 
condition is very important for simulations with WDMB. 
Star formation occurs then on a time scale 

^ _ tdyn _ 1 1 

governed by the dynamical time scale tdyn and a star forma- 
tion efficiency (SFE) e. The SFE is quite small in the Galac- 
tic disc at present, at most a few percent (e.g. Silk 1997). 
Therefore, most simulations are run with a SFE e — 0.01. 
But to account for the possibility of more efficient star for- 
mation in the early universe, we also use e = 0.1 and e = 1 
for some of the simulations. As is customary, the actual point 
in time, when a SPH particle is converted to a star parti- 
cle, is then determined probabilistically, based upon its star 
formation time scale. 

The SFE is kept constant in time, which is of course 
unrealistic for the cases with high SFE e — 0.1 and e = 1. 
The high-SFE simulations were primarily run to see what 
effect the SFE has on the onset of star formation at high 
redshifts. 

Feedback is not included in the simulations. Even 
though it may be important in the early stages of 
galaxy formation (e.g. Sommer-Larsen, Gotz fc Porti- 
nari 2003), it is probably not at later times. In the 
warm dark matter (WDM) cosmological scenario, which 
is used here (see below), the angular momentum prob- 
lem, as well as other problems on galactic scales re- 
lated to CDM, c an be significantly reduced wit hout in- 
voking feedback dSommer-Larsen fc Doleovl l200ll) . In any 
case, the absence of feedback in the simulations can be 
regarded as a conservative approach, since its inclusion 
would possibly result in even more real istic galactic discs 
iSommer-Larsen. Gotz fc Portinarill2003|) . 



5.3 Cosmological Parameters and Initial 
Conditions 

The cosmological initial conditions are based on a flat 
AWDM model with Qm = 0.3 and = 0.7. The Hub- 
ble parameter is chosen to be Ho ~ lOO/i km / s / Mpc = 
65 km / s / Mpc, yielding a present age of the universe of 
14.5 Gyr. The initial power spectrum of density fluctua- 
tions is cluster abundance normalized to give a erg — 1.0 
at z=0 (Eke, Cole fc Frenk 1996 — adopting ag = 0.9 
does not lead to any signiflcant changes of the results pre- 
sented in the following). For the WDM particle, we use a 
free-streaming mass of 1.5 • 10^^ Mq, corresponding to 
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Table 2. Physical properties of the two selected disc galaxies in the different simulation runs. A '0' in the column '# of clumps' means 
that the simulation was run with evenly distributed baryonic matter, otherwise the number of clumps in the initial galaxy forming region 
is given. [Fe/H] is the metallicity, with 'prim.' denoting the primordial value, e the star formation efficiency, M« and Meg the masses in 
stars and cold gas (T < 3 • 10* K), respectively, Vc the circular velocity, j'tot the specific angular momentum of the disc, including stars 
and cold gas, and j» the specific angular momentum of the stars only. The quantities in the last five columns are for within 40kpc of 
the center of the disc. 



Galaxy 


# of 


[Fe/H] 


£ 




Meg 


Vc 


itot 






clumps 






[lOi^Mo] 




[km 1 s] 


[kpc ■ km / s] 


[kpc ■ km / s] 


19 





prim. 


0.01 


4.97 


0.83 


209 


511 


341 


19 





-0.5 


0.01 


6.98 


1.81 


218 


1433 


961 


19 


248 


prim. 


0.01 


5.41 


0.76 


213 


557 


422 


19 


248 


prim. 


0.1 


4.64 


0.05 


195 


271 


270 


19 


248 


prim. 


1.0 


4.35 


0.004 


185 


292 


292 


19 


248 


-0.5 


0.01 


7.98 


1.27 


229 


1119 


878 


27 





prim. 


0.01 


3.50 


0.53 


179 


466 


224 


27 





-0.5 


0.01 


5.17 


1.01 


196 


934 


493 


27 


177 


prim. 


0.01 


3.65 


0.62 


182 


536 


331 


27 


177 


prim. 


0.1 


3.70 


0.27 


178 


458 


336 


27 


177 


prim. 


1.0 


3.34 


0.09 


161 


432 


404 


27 


177 


-0.5 


0.01 


6.31 


1.13 


217 


851 


621 



a free-streaming scale of O.llft ^ Mpc and a WDM particle 
mass of \.2]V''' ^ keV. These numbers are consiste nt with the 
values used by ISommer-Larsen fc Doleovl ll200lll for solving 
(to within a factor of two) the angular momentum problem 
with WDM. In addition, the same sharp cut-off of the WDM 
power spectrum below the free-streaming scale is used here. 

In order to select interesting candidate galaxies for fur- 
ther study, we ran at first a dark matter-only N-body sim- 
ulation in a lOfe"^ Mpc box with 128^ particles , using the 
Hydra code iCouchman. Thomas fc Pearcelll^95^ and start- 
ing at an initial redshift Zini = 19. The usual set-up, where 
particles are displaced from a regularly spaced grid, can not 
be used in WDM simulations, as this leads to the formation 
of spurious low-mass haloes evenly spaced along the form- 
ing filaments, which are non-physical artifacts. (See Gotz & 
Sommer-Larsen 2002 for a detailed disc ussion.) We therefore 
start from glass-like initial conditions JWhitd ll9961 where 
the particles are irregularly distributed, but still (almost) 
evenly spaced. 

At redshift z = 0, haloes in the pure dark matter simula- 
tion are identified with a friends-of-friends algorithm, where 
gravitationally unbound particles are removed. The galax- 
ies then selected for detailed simulation with the TreeSPH 
code had to be sufficiently isolated (at least 1 Mpc away 
from galaxy groups and 0.5 Mpc away from larger galaxies 
at z = 0) and had not to undergo a major merger (with 
mass ratio more than 1:3) since z = 1. 

To simulate the formation of the selected galaxies, bary- 
onic matter has to be added to the initial conditions. We 
aimed at obtaining a baryonic fraction /b ~ 0.1, consistent 
with nucleosynthesis constraints and with o bserved baryonic 
fractio ns in groups and clusters (see e.g. lEttori fc Fabiarj 
il999l) '). Our choice of /b is a bit low compared to recent 
estimates of and VLm (section 4), but this has no conse- 
quences for any of the following results. 

For the simulations with baryonic clumps, we choose 
the most extreme situation where all of the baryonic matter 



is initially contained inside spherical clumps, whereas the 
dark matter is exclusively found outside of them. That way, 
we maximize the effect, which the clumps have on galaxy 
formation. The easiest way to achieve the desired baryonic 
fraction /b ~ 0.1 is then to remove 10% of the original 128'' 
dark matter particles from spherical subregions of the cos- 
mological simulation, and fill these regions with gas parti- 
cles. To get (approximately) the same number of particles as 
in the simulations, where the baryonic matter is distributed 
evenly, we use 128"' SPH particles for that purpose, such 
that each has a mass /b times that of a dark matter parti- 
cle. We therefore get in the clumps a lOx higher resolution 
(in terms of particle number density) than in the surround- 
ings, where only dark matter is found, and can thus study 
the interaction of the clumps in great detail. 

Because we are constrained by the fact, that there has 
to be an integer number of clumps, each containing an in- 
teger number of SPH particles, we ended up with placing 
128^^ — 16384 baryonic clumps into the box, each with 128 
baryonic particles and replacing 13 dark matter particles. 
This gives a baryonic fraction /b = 0.102, very close to the 
desired value. The initial conditions with the clumps contain 
then (1 - /b) • 128^ = 1884160 dark matter particles (with a 
mass of niDM = 4.0- lO'^/i'^M© for the 10/i~^ Mpc box) and 
128^ SPH particles (with a mass of mb = 4.0 ■ IffhT^Mo). 
The mass and radius of a clump are Ma = 5.2 • lO**fe"^M0 
and i?ci = O.llft"'^ Mpc, respectively, with an average sep- 
aration of del = 0.39/1"^^ Mpc (all length scales given are in 
comoving coordinates). Our choice of clump mass is moti- 
vated by the desire to have Md <^ Mfs, but still large enough 
to be resolved in the simulations. Any other choice of Md 
would lead to very similar results, provided Md <^ Mfs — 
see below. 

The clumps are distributed irregularly throughout the 
simulation volume, with their centers forming a glass-like 
structure. The 13 closest dark matter particles to a clump 
center are removed, and the resulting void is filled with a 
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sphere containing 128 SPH particles, also with a glass-like 
distribution, which furthermore is oriented at random. This 
procedure assures that the initial density field is not dis- 
turbed, and that neighboring clumps do not overlap initially. 
(Random placement of the clump centers would result in 
many overlapping clumps as their average separation is not 
much larger than their size, dd ~ 3.5-Rci.) 

For the control simulations with evenly distributed 
baryonic matter, one SPH per dark matter par- 
ticle was added within one SPH softening length 
iSommer-Larsen fc Dolgovll200j) and with a mass /b times 
the original dark matter particle mass. Since no dark mat- 
ter particles were removed, their masses were reduced to 
(1 — /b) the original value. Thus rrih is the same as in 
the simulations with clumps, whereas now niuM ~ 3.6 ■ 
10'' h-'^MQ. 

From these initial conditions, which now contain 
SPH particles, regions were cut out corresponding to the 
subvolumes, from which the galaxies of interest form in the 
pure dark matter simulation, as identified by tracing their 
particles back to the initial conditions. (Care is taken to 
make sure that no clumps are intersected by the region 
boundaries in simulations with baryonic clumps.) Only 
inside these areas are the gas particles included, whereas for 
the surrounding region the original pure dark matter initial 
conditions were used, with the dark matter particles in- 
creasingly coarsely resampled with increasing distance from 
the galaxy forming subvolume. That way, the formation of 
individual galaxies can be followe d in detail within the cor - 
rect cosmological setting, se e e.g. lNayarro fc Whit 3 il994): 
iGelato fc Sommer-Larsen| l| l999|): [ T hacke ^^^ouchmanl 
j200d) : ISommer-Larsen. Gotz fc Porti nari (20qJ! 

Initially, the SPH particles are assigned a thermal en- 
ergy corresponding to a temperature T ~ 55 K. Gravita- 
tional interactions are so f tened according to the prescrip- 
tion of iHernauist fc Katd il989h . The softening lengths are 
1.3h~^ kpc for gas and star particles, 2.9h~^ kpc for the dark 
matter particles inside the galaxy forming region in simula- 
tions with clumps, and 2.8h~^ kpc for such dark matter par- 
ticles in runs with evenly distributed baryonic matter. The 
softening lengths of the dark matter particles in the coarsely 
resampled surrounding are correspondingly larger. 

5.4 The Simulations 

We selected two dark matter haloes (numbers 19 and 27) 
from the cosmological N-body simulation for further study 
with the TreeSPH code. The haloes have characteristic cir- 
cular velocities comparable to or a little below that of the 
Milky Way. These galaxies were chosen according to the cri- 
teria laid out above, and because they show distinct disc- 
like morphologies and exponential surface density profiles 
in the TreeSPH-runs with evenly distributed baryonic mat- 
ter. For each of the haloes, six simulations with different 
combinations of cooling function metallicity, star formation 
efficiency, and presence or absence of baryonic clumps, were 
run, as can be seen in Table 15.31 The simulations of the 
larger halo, #19, contain about 34000 dark matter and 
22000 (without clumps) or 32000 (with clumps) SPH parti- 
cles initially, whereas the runs for halo #27 consist of about 
26000 dark matter and 14000 (without clumps) or 23000 
(with clumps) gas particles at the start. The significantly 
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Figure 7. Surface density profiles for Milky Way like disc galaxy 
#19 in stars (top) and cold gas with T < 3 ■ 10* K (bottom) at 
2 = 0. Thick lines denote runs with evenly distributed baryonic 
matter, and thin lines those with baryonic clumps. Simulations 
with primordial abundance are shown as solid lines, and the 1/3 
solar metallicity cases ([Fe/H] = —0.5) by dashed lines. The SFE 
is £ = 0.01, unless otherwise noted. The scale lengths for the 
stellar discs lie between 1 kpc and 1.5 kpc. 

higher number of SPH particles in simulations with bary- 
onic clumps is due to the fact that their initial regions have 
to be extended to avoid that clumps are intersected by the 
boundaries of the galaxy-forming subvolumes. 

5.5 Results 

5.5.1 Properties at z — 

In all of the TreeSPH-simulations, the galaxies show clear 
disc-like morphologies and kinematics, with the bulk of the 
stars on approximately circular orbits in the disc. Most of 
the remaining stars are found in an inner, bulge-like compo- 
nent and finally a small fraction in a dynamically insignif- 
icant stellar halo surrounding the disc. The disc galaxies 
formed in our simulations are thus qualitatively very similar 
to observed disc galaxies, like the Milky Way, irrespective of 
the presence or absence of baryonic clumps initially, and of 
the assumed values for the IGM metallicity and SFE. 

All of the discs have approximately exponential surface 
density profiles in stars and cold gas out to four scale lengths. 
As an example, the profiles for the different simulations 
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Figure 8. The normalized specific angular momenta in stars 



at 



of the disc galaxies formed in the simulations. 



Small symbols are for runs with a SFE e = 0.01, medium sized 
ones have e = 0.1, and the large ones e = 1. The round symbols 
denote primordial abundance, and the triangles a metallicity of 
[Fe/H] = —0.5. Open symbols are used for the control simulations 
with evenly distributed baryonic matter, whereas filled symbols 
are for runs with baryonic clumps. The solid line s hows the me - 
dian value from observati onal data on disc galaxies Byuri|| l992l) , 
obtained as explained in ISommer-Larsen fc Dolgov J200lft . but 
using Ho = 65 km / s / Mpc, the dashed and dotted lines bracket 
the Icr and 2cr intervals around the median. 
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of galaxy #19 are shown in Fig. The exponential scale 
lengths for the stellar discs he between Ikpc and 1.5 kpc. 
The runs with 1/3 solar metallicity (dashed lines in Fig. |7|l 
show somewhat flatter gas and stellar profiles and a lower 
density in cold gas, since the more effective cooling leads to 
a higher conversion rate of gas into stars. The simulations 
with a higher SFE (as marked in Fig.UJ have, as one would 
expect, much depleted gas discs, and lower surface densities 
in stars, since their star formation rates peak at higher red- 
shifts (see below). Thus more stars are formed initially in 
the halo, and less gas is left to form the disc stars. 

The scale lengths of the stellar discs (see Fig. |7|l are 
somewhat low compared to observed values — thus the an- 
gular momentum problem has not been completely over- 
come. This can also be seen by looking at the normal- 
ized specific angular momenta in stars J« = j*/V^ plot- 
ted in Fig. IHl for the two discs in all the different simu- 
lation runs. As argued bv ISommer-Larsen. Gelato fc VedeJ 
il999t) one expects J* to be almost independent of Vc 
on both theoretical and observational grounds. On av- 
erage, the normalized specific angular momenta of our 
discs lie only about a factor of two below the median 
observ ed value, as obtained by ISommer-Larsen fc Doleovl 
l|200ll) . This compares well with other attempts to solve 
the angular-momentum problem , e.g. the WDM models of 
ISommer -Lar^^n ^ D^l^yyl l|2(^(^l [l and CDM scenar ios with 
feedback~ ilSoimnB^Lar'seirG5rz fc Portinaril l2003l) , which 



Figure 9. Star formation rates in Mq/ yr for the different sim- 
ulations of galaxy #19, with a magnification of the high-redshift 
tail at the bottom. Line styles are as in Fig. 171 

also come within a factor of two to three of the observed 
specific angular momenta. 

In Table 15.31 and Figs. Q and |5] we have included the 
results for the simulations with high SFE e = 0.1 and e = 
1 at = 0, even though these runs were initially carried 
out only to study the efl'ect, a possible high SFE at large 
redshifts would have on early star formation. As expected, 
these simulations produce discs which have an unrealistically 
low cold gas fraction, indicating that the SFE can not stay 
constant at these large values. In spite of that, their specific 
angular momenta lie in the same range as that of the other 
runs (see Fig. ISJ, indicating that SFE is not an important 
parameter in WDM models with baryonic clumps for solving 
the angular momentum problem. 

5.5.2 The High-Redshift Evolution 

The present-day properties of disc galaxies are thus unaf- 
fected by the presence or absence of baryonic clumps. In 
both cases do we find discs at redshift z — with the same 
kinds of profiles and specific angular momenta. The reason 
for this lies in the behavior of the clumps at the beginning 
of structure formation. 
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Figure 10. Same as Fig. |2l but now for galaxy #27. 



higher density. If one only applied a density criterion as the 
star formation recipe, one could get star formation started 
very early in WD MB. But the initial gravitational compres- 
sion and deformation of the clumps heats the baryonic ma- 
terial, and afterwards a balance between Hubble expansion 
and tidal gravitational squeezing is achieved for some time, 
keeping the baryonic density in the clumps approximately 
constant. Only after turn-around has the baryonic material 
in the clumps a chance to cool efficiently and to start form- 
ing stars. This is why we have included the condition, that 
the gas velocity field needs to be compressive in star forming 
regions, into our star formation recipe. 

The first stars form at an epoch 2,4 which is of or- 
der a period t* after the first region with a compressive 
gas velocity field and nu > n-H,c appears. For the standard 
WDM simulations with e — 0.01 z*,! ~ 4.5; cf. Figures 9 
and 10 (the [Fe/H]=-0.5 cooling function simulations are 
clearly inappropriate when discussing the formation of the 
first stars). The similar WDMB simulations have ~ 6.5. 
The WDMB simulations with higher SFEs (e = 0.1 and 1) 
have 2*4=7-9.5. We did not run standard WDM simulations 
with e = 0.1 and 1, but based on the e — 0.01 simulations we 
estimate that even for such high SFEs 2*4 < 5.5. Hence going 
to the WDMB model brings the warm dark matter model in 
much better agreement with observations of high-2 galaxies 
(e.g. Hu et al. 2002 — 2=6.56) and QSOs (e.g. Fan et al. 2001 
— several with 2 > 6). In fact for the high SFE simulations 
there is even marginal agreement with the value of the red- 
shift of reionization o f the Universe z^e > 10 recently claimed 
by the WM AP team dSoergel et al.l2o63() , provided reioniza- 
tion happens shortly afte r the onset of star formation — but 
see be low. We note that ISommer-Larsen. G5tz fc Portinaril 
i2003l) find that the SFE likely was considerably higher dur- 
ing the early phases of galaxy formation than at present. 
For 2 > 10 the gas velocity field is everywhere expansive, so 
it will not be possible to push 2,4 beyond «10. 



When the first sheets and filaments form at high red- 
shifts, the baryonic clumps get rapidly distorted into el- 
lipsoidal shapes by the tidal field of these structures. The 
clumps get 'gently' incorporated into the filaments, and 
there are no shocks occurring during this, or when the 
clumps eventually physically touch. After 1.0-1.5 Gyr {z—4- 
5), the filaments in the simulations with baryonic clumps 
are almost evenly filled with baryonic matter, after which 
the same processes of galaxy formation take place as in 
the ordinary WDM scenario with evenly distributed bary- 
onic matter (The clumps do though survive much longer in 
the voids, but this is not of interest here, since the early 
galaxy/star formation in the simulations takes place in the 
filaments). During filament formation, also the dark matter 
becomes distributed smoothly in the filaments. This general 
behaviour strongly suggests that the main results of the sim- 
ulations do not depend of the specific choice of clump mass, 
as previously mentioned. 

But there is a significant difference between the clump 
and no-clump models in terms of the early star formation 
history. Figs. |^ and 1101 show the star formation rates for 
the two selected galaxies in all the model runs. In the sim- 
ulations with baryonic clumps, the first stars form earlier 
than in those without them. This is not surprising, since the 
baryonic matter in the models with clumps starts out with a 



6 SUMMARY AND CONCLUSIONS 

We have investigated the implications of strongly inhomoge- 
nous, primordial baryon distribution on sub-galactic scales 
{WMq < Af < 10" Af©) for Big Bang Nucleosynthesis, 
CMB anisotropies and galaxy formation in the context of 
the warm dark matter model. Big Bang Nucleosynthesis is 
only slightly changed relative to SBBN, but the change in 
recombination history at 2 ~ 1500 — 700 relative to "stan- 
dard" theory leads to differences in the anisotropy and po- 
larization power spectra, which should be detectable by the 
Planck satellite provided systematic efi'ects can be under- 
stood. We show by fully cosmological, hydro/gravity sim- 
ulations that the formation of galactic discs is only weakly 
affected by going from the WDM to the WDMB scenario. In 
particular, the final disc angular momenta (at 2=0) are as 
large as for the standard case and the "disc angular momen- 
tum problem" is solved to within a factor of two or better 
without invoking (hypothetical) energetic feedback events. 
A very desirable difference from the standard WDM model, 
however, is that the on-set of star (and AGN) formation 
happens earlier. For the "optimal" free-streaming mass of 
Mf ~ 1.5 • 10" /i"^ Af© the redshift of formation of the 
first stars increases from 2,,i=4-5 to ^6.5, in much better 



© 2003 RAS, MNRAS OOO.nTOl 



High-z Galaxy Formation with WDM 13 



agreement with observational data on high-redshift galaxies 
and QSOs. The WDMB model is, however, for this choice of 
free-streaming mass still limited to < 10, so probing the 
"dark ages" with JWST, ALMA etc. wiU offer direct tests 
of this theory. 
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